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ABSTRACT 


A fundamental study of the stability of a dynamic 
gas-liquid interface between rotating cylinders is reported. 
The study was initiated for the purpose of seeking factors 
which have a significant role in the process of gas inges- 
tion, or gas entrainment, in viscoseals . 

The simplified model of smooth, cylindrical surfaces 
was selected for mathematical tractability and to provide 
a visual study, using a transparent acrylic housing, without 
the obscurity of the more complex fluid flow resulting from 
the presence of the grooved surfaces employed in viscoseals . 
The visual study was supplemented by employing stroboscopic 
photography and high-speed motion picture photography. 

A phenomenological mechanism of gas ingestion was 
established, theoretically and experimentally. It was found 
that gas entrainment can result from a gas -liquid interface 
instability caused by a velocity of a portion of the inter- 
face toward the more viscous fluid and/or an acceleration 
of a portion of the interface toward the more dense fluid. 
Results of the study indicate that surface tension tends to 
stabilize the interface and prevent or delay gas ingestion. 

It is concluded that gas ingestion in viscoseals is 
an essentially inherent phenomenon which cannot be readily 
eliminated through optimization of the viscoseal geometry. 
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A brief discussion is given on the implications of 
the interface stability study in the operation of journal 
bearings, face seals, and buffered bushing seals, in 
addition to a discussion of its application to viscoseals. 

Eccentricity, shaft runout, and vibrations are 
credited with contributing to conditions necessary for the 
instability of an interface in rotating machinery . 
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CHAPTER I 


INTRODUCTION 

A fundamental study of the stability of a dynamic 
gas -liquid interface was initiated for the purpose of 
providing insight into the mechanism of gas ingestion in 
viscoseals . 

A viscoseal is a dynamic shaft seal which develops 
a pressure gradient in the axial direction within the 
annular space between a shaft and housing having a relative 
angular velocity and normally not in contact with each 
other. Helical grooves are provided on the shaft and/or 
the housing, and the relative motion of the two surfaces 
causes a viscous pumping action on the fluid in the annular 
space, developing the desired pressure gradient in the 
axial direction and establishing a gas -liquid interface 
within the annulus. A typical viscoseal, with a grooved 
shaft and smooth housing, is shown in Figure 1. 

Since the viscoseal does not depend on rubbing con- 
tact in providing sealing action, it offers the possibility 
of low wear rates and essentially zero-leakage rates for 
extended periods of operation. These considerations are 
important in space applications and nuclear operations, and 
also make the viscoseal an attractive contender for many 
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Grooved Interface 

Shaft 


Figure 1. Elements of a typical viscoseal. 
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other industrial applications , particularly those involving 
high-speed rotating machinery. 

The viscoseal is essentially the same physical 
apparatus as a screw extruder, a device of earlier origin 
operated for the primary purpose of producing a flow of a 
very viscous liquid or a plastic, usually under pressure. 

The first patent of an extruder employing a screw was taken 
out in 1879 by an Englishman named Gray [ 1 ] ^ . John Royle 
developed a screw machine in the U.S. at the same time. 

Lawaczeck [2] conducted tests with the viscopump in 
1916 in an extension of the screw extruder concept to the 
performance of the viscoseal function of generating an 
axial pressure gradient while producing no net flow through 
the annulus . 

An experimental screw oil pump, reported by Pearsall 
[3] in 1924, used a threaded shaft to produce a high-pressure 
oil flow, and also employed screw threads "to act as a gland 
and prevent the oil escaping" along a sleeve. Pearsall 
derived a relation between the volumetric flow rate and the 
outlet pressure. He commented that "when using such a 
thread to retain oil and prevent it leaking out the end of 
the bearing, it has been found in practice necessary to run 


■^Numbers in brackets refer to similarly numbered 
references in the bibliography. 
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at a speed about two and a half times that given by the 
equations for no flow." 

Rowell and Finlayson [4] presented an analysis of 
the screw viscosity pump in 1928 and commented on its ’‘most 
immediately interesting application" in the field of bearing 
seals . 

Interest in the viscoseal concept continued, but the 
subject received only a limited amount of attention until 
the needs for advanced sealing concepts in the space and 
nuclear programs gave impetus' to the development of the 
viscoseal at a more rapid pace. 

In 1963, McGrew and McHugh [5] reported an observed 
phenomenon which they called "seal breakdown" which appeared 
as a leakage of about one drop of liquid every three minutes 
and a decrease in the seal pressure. In their test rig 
arrangement, the axis of the shaft was parallel to the 
direction of the force of gravity. The liquid was above 
the interface and the leakage which was encountered was in 
the form of liquid dripping out the open end of the housing 
to the atmosphere. They studied the relation between the 
beginning of seal breakdown and the Reynolds number, Froude 
number, cavitation number, and Weber number. The Weber 
number, involving the fluid properties of density and sur- 
face tension, appeared to provide the best correlation with 
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seal breakdown, but was not conclusive. Three silicone 
fluids and water were tested. Seal breakdown for three of 
these fluids occurred in a range of Weber numbers from 
800-1100. For the remaining fluid, which was one of the 
silicone fluids, the Weber number was 220 at the beginning 
of seal breakdown. Seal breakdown occurred under laminar 
flow conditions for another of the silicone fluids, while 
breakdown came in the transition range for one of the fluids 
and in the turbulent range for the other two fluids. 

In a later experimental study, King [6] also employed 
a vertically mounted shaft with the liquid above the inter- 
face and observed a "secondary leakage brought on by splash- 
ing and instabilities at the liquid-vapor interface." Liquids 
tested by King were oil, water, and potassium. In the tests 
with water, the point at which the secondary leak could first 
be detected varied with the speed and radial clearance. At 
0.0025 in. radial clearance, it was first noticed at 12,000 
to 14,000 rpm; at 0.0029 in., it was noticed at 7,000 to 
8,000 rpm; at 0.0044 in., it was noticed at the starting 
speed of 4,000 rpm. The Weber numbers corresponding to 
these points were 1290 , 5 72 , and 214, respectively. With 
a 0.0015 in, clearance, the shaft was operated at a speed 
up to 14,000 rpm with a corresponding Weber number of 883, 
but no leakage was noticed. 
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In experimental work reported by Stair [7] in 1965, 
the viscoseal shaft was oriented horizontally. During the 
course of this experimental work no evidence of a seal 
break was observed, but a phenomenon called air ingestion 
occurred. When the shaft speed was increased to a suffi- 
ciently high value, air bubbles were observed to rise in 
the transparent pressure tap lines, and bleeding the air 
from the lines would not stop the ingestion. The axial 
pressure gradient decreased and caused the air-water 
interface to move . toward the atmospheric end of the seal, 
but there was no observed leakage except when the interface 
moved outside the annular space. At increased shaft speeds 
the amount of air ingestion would increase and in some 
cases the shaft reached a speed at which the pressure profile 
would become unstable and severe pressure pulsations were 
observed. Stair concluded from his data that air ingestion 
increases with shaft eccentricity and is reduced by using 
wider grooves. He noted [8] that no air ingestion was 
detectable with laminar operation of the seal . 

An experimental study of gas ingestion was reported 
by Ludwig, Strom, and Allen [9] in March 1966. Test fluids 
used were water and sodium. In the tests with water, a 
transparent housing was used to permit visual observations 
of the gas ingestion region. For each fluid, tests were 
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made using a grooved shaft, and using a grooved housing, 
with the same groove depth, width, etc. used for each. Gas 
ingestion was found to be more pronounced with the grooved 
shaft, and the authors attributed the difference to the 
centrifuge action forcing bubbles into the grooves where 
they are pumped to the high-pressure end of the seal. 

Although increasing attention has been given to the 
problem of gas ingestion, the basic mechanism causing it 
is not well understood. It is not understood how gas 
bubbles which become mixed with the liquid are transported 
to the high-pressure end of the seal. A separate question 
concerns the mechanism by which the interface becomes 
interrupted to permit the gas bubbles to become entrained 
in the liquid. 

While the mechanism of seal breakdown or secondary 
leakage may be significantly different from that causing 
gas ingestion, both appear to result from some form of 
instability of the gas- liquid interface, suggesting the 
need for the study undertaken here. 

A direct study of the stability of the dynamic 
interface actually encountered in the viscoseal presents 
an extremely difficult problem. Therefore, for the purpose 
of this investigation, a simplified model leading to a more 
basic study of the stability of a dynamic interface was 
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selected as a matter of practical expedience. The model 
which was chosen as having a significant similarity to the 
flow in a viscoseal, but sufficiently simple to permit an 
analytical study, was a gas and liquid confined in the 
annular space between the smooth surfaces of a stationary 
housing and a rotating cylinder. This model is the same 
as the viscoseal shown in Figure 1, page 2, except that 
a shaft with a smooth surface is substituted for the grooved 
shaft , 



CHAPTER II 


ANALYTICAL STUDY 
I. MATHEMATICAL MODEL 

The mathematical model for the analytical study was 
based on further simplifications of the flow model. Because 
of the small clearance between the surfaces having relative 
motion, the effects of curvature were considered to be of 
small magnitude in comparison with other variables of 
interest. Therefore, the model for mathematical analysis 
was obtained by assuming an interface between two incom- 
pressible, Newtonian fluids confined between two infinite 
flat plates moving at equal but oppositely directed 
velocities U/2 in the x-direction , as shown in Figure 2. 

The theoretical analysis was approached by consider- 
ing separately the effects of the fluid viscosities and the 
effects of fluid acceleration. The model for the analysis 
treating viscous effects was patterned after the work of 
Saffman and Taylor [10] in a study of the stability of an 
interface separating two viscous fluids confined between 
two closely spaced parallel plates with the fluids having 
a velocity perpendicular to the ' interface . The plates were 
assumed to be stationary and inertia effects were neglected 
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Figure 2. Model for analysis of interface separating 
two fluids confined between parallel plates having relative 
motion parallel to the interface. 
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in this analysis. The analysis taking acceleration effects 
into account was developed along lines similar to the work 
of Taylor [11] in a stability analysis of an interface 
separating two fluids of different densities when the fluids 
are accelerated in a direction normal to their mutual inter- 
face.. His analysis was not restricted to a narrow space 
between two plates, and viscous effects were not included 
in the analysis. 

II. STABILITY ANALYSIS - VISCOUS EFFECTS 


For the Cartesian coordinate system shown in Figure 


2, the Navier-Stokes equations are 




y v 2 u; 


CD 


p 


Dv _ 
Ft " 



+ yv v; 


( 2 ) 


P 


Dw 

Dt 



2 

y V w. 


(3) 


The clearance c between the plates is assumed to be 
small and the fluid motion is assumed to be everywhere 
parallel to the plates. Derivatives of velocity components 
in directions other than normal to the plates are assumed 
small in comparison with derivatives along the normal, and 
are neglected. If the acceleration terms are assumed to be 
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small, the equations of motion become 


and 


0 = 


3P , 3 2 u 

** * 77 ’ 


(4) 


0 = - 


( d? \ 3^V 

W + pg ) + 77 


(5) 


Integrating equations (4) and (5) with respect to z, 


3u = HP z + c 

3 Z li 3X Z 1’ 


( 6 ) 


and 


3 v 
3 z 


■ Hly * pg ) 2 + c ; 


(7) 


A second integration with respect to z yields 


u 


3P_ 2 

2v 3 x z 


+ C-^z + Cj , 


and 


( 8 ) 


v 


27 (If + pg ) z2 + C 2 Z + C 4 ' 


(9) 


The appropriate boundary conditions are, 
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at z = 0: u=-y; v = 0 ; 

( 10 ) 

at z = c: u = + ^; v = 0 „ 

Substituting the boundary conditions into equations (8) and 
(9) , the constants are 


C4 0 , 

C = u - J, 21 c 
L 1 c 2u 3x c ’ 


( 11 ) 

( 12 ) 

(13) 


C 


2 


1 

27 




c . 


(14) 


Equations (8) and (9) become 


_1 3P _2 . Uz c 3P _ U 
■ u 


u 2u 3x z + c 2 m 3x z " 2* 


and 


v 


1 /3P A \ 

" 27 V3y pg J z 


2 c /3P . \ „ 

- 27 1 3y + pg J z 


The mean velocities are 



u - — / u dz , 

c 7 7 


v = i / v dz , 
c /O 


(15) 


(16) 


(17) 

(18) 
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u = — 


v = 


u = 


V = 


1 

C 


1 

3P z 3 

+ Uz2 

c 3P z 2 Uz 

C 

(19) 

2u 

9x 3 

2c 

2y 9x 2 2 

9 

0 

i 

2y 

fli + 

\9y 

x 3 
P gj ^ — 

2 “I c 

c rap . \ z z 

" 2TT lay + r-J 0 > 

(20) 

c 2 

12y 

3P 

3x ’ 




(21) 

c 2 

12y 

fi L + 

\3y 

Pg) • 



(22) 


The form of the equations for the mean velocities 
permits the introduction of a velocity potential <j> , given 
by 

2 


u = - 


9P 

12 y 3x 


3(f> 

3x » 


(23) 


and 


c 2 ( 3P x \ _ 3(f) 
“ — ' ~ + PgJ = “ 


12 u V 9y 


3y 


(2 4) 


The interface is assumed to have a slight wave -like 
disturbance of wavelength 

L = 2tt/K, (25) 

described by 

S = Be mt cos Kx, (26) 

which is the equation for the surface of separation of the 
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two fluids, and which gives the displacement of the inter- 
face from the mean position of the interface, as shown in 
Figure 3. 

The stability problem becomes one of determining the 
algebraic sign of m in equation (26) . If m is positive, 
the magnitude of the interface disturbance grows exponen- 
tially with time and the interface is unstable. Neutral 
stability exists when m is zero, and the interface is stable 
to a small disturbance if the sign of m is negative. 

On the interface, the y-components of the velocities 
of the two fluids are equal, and are given by 


3 *! 

W 


9 4 2 

9y 


— 9S 

v + at 


(27) 


= v + mBe mt cos Kx. (28) 

The solution of the continuity equation, 

V 2 <j> = 0, (29) 

must satisfy equation (28) , and the disturbance velocity 
must vanish at infinity. The appropriate solutions are 

+1 - vy - 2® e mt - K r cos Kx, (30) 

and 


• _ — . mB mt + Ky 

<f >2 vy + — £ e 7 


cos Kx. 


(31) 



i 


Fluid 1 


cos Kx 


Fluid 2 


Figure 3. Model for analysis of a disturbed 
interface separating two viscous fluids moving transverse 
to the interface and confined between parallel plates 
having relative motion parallel to the interface. 






17 


Integrating equation (24) with respect to y, and 
substituting the velocity potentials from equations (30) 
and (31) , the solutions for the pressures are 

12 

P, « - —^(vy - e rat ' Ky cos Kx) - P lg y + (v , (32) 

c 

and 

P 2 = " — e mt + Ry cos Kx) - P 2 gy + . (33) 

c 

The difference between the pressures on the two 
sides of the interface will depend on the curvature of 
the interface [12], and may be expressed by the Laplace 
equation , 

( p i - v - t * y > (34) 

where T is the surface tension and R-^ and R 2 are the radii 
of curvature of the interface in mutually perpendicular 
planes. The interface curvature may depend in a rather 
complicated manner on such factors as the distance sepa- 
rating the plates, surface roughness, the materials of the 
solid surfaces, the liquid, and interface dynamics. For 
example, the contact angle for a receding interface would 
be expected to be less than the angle of contact of an 
advancing interface [13] . However, it is assumed here for 
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simplicity that the angles of contact with the solid 
surfaces are zero and that the radius R^ is half the 
distance separating the solid surfaces, or 


1 = 1 
R x c/2 


( 35 ) 


Because the interface disturbance was assumed to be small, 
the curvature in a plane parallel to the plates may be 
approximated by 



(36) 


When equations (35) and (36) are substituted inti 
equation (34) , the pressure difference becomes 


(P, - P 9 ) = T 




(37) 


Substituting equation (26) into equation (37) and setting 
equation (3 7) equal to the pressure difference across the 
interface from equations (32) and (33) gives 


Tf 2 *,^2 mt „„ _ 12vS r ^ 

T(— - BK e cos Kx) = - — -^--(p^ ” v 2' 

c 


12m B r „mt - KS A „mt + 

+ — s— (p,e + p 9 e ) cos Kx 

c 4 K 1 z 


- CPj - P 2 )gS + (C s - C 6 ) 


(38) 
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When S, as defined by equation (26), is zero, 


equation (38) reduces to 

T(|) = (C 5 - C 6 ). (39) 

Substituting this into equation (38) , then dividing through 
by S and multiplying through by K = 2 tt/L gives 


12m, Q 


-KS + KS-. 2 it 

+ y e ) = r - 


12 v 


Cu ^ u 2 ) + (p 1 ^2^® 


8tt 3 T 


(40) 


A disturbance is unstable when the right-hand side 
of equation (40) is positive, since the left-hand side is 
positive when m is positive. Therefore, the first term on 
the right-hand side of the equation must be positive and 
numerically larger than the second term on the right in 
order for a disturbance to be unstable. 

Thus, surface tension tends to cause the interface 
to be stable to a small disturbance, and an interface 
velocity toward the more viscous fluid contributes toward 
the development of an instability. 

Information concerning the wavelength L u of unstable 
disturbances can also be obtained from equation (40) . Set- 
ting the right-hand side of equation (40) to be greater 
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than zero, 

L u 2 p| Cu x - y 2 ) + (P X - P 2 )gJ - 4 tt 2 T > 0, (41) 

which is the case for unstable disturbances, and solving 
for L u gives 


L u > 


2ttcT 


1/2 


Jl2v(i 


y?) + c (p, - p 9 )g 


] 


1/2 


(42) 


When the necessary conditions exist, as given by equation 
(40), for small disturbances to be unstable, the wavelength 
of such unstable disturbances must be greater than the 
critical value given by equation (42) . This value is pro- 
portional to the square root of the surface tension, so for 
the hypothetical case in which the surface tension is zero, 
the wavelength of unstable disturbances would not be limited 
to a minimum value. The result is the same as that found 
by Saffman and Taylor [10], even though the plates were 
considered to be stationary in their analysis. 

That an unstable disturbance will continue to grow 
when it reaches a finite magnitude can be shown readily 
by a qualitative determination of the relative magnitude of 
pressure gradients at specific locations in Figure 4, where 
constant pressure lines in the liquid are sketched, but not 
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necessarily to scale. For the purpose of the qualitative 
analysis under consideration, assume that the gravitational 
force is normal to the plane of the plates in Figure 4. 

Since the primary concern here is with a gas -liquid inter- 
face, and liquid viscosities are usually much larger than 
the viscosity of gases, the viscosity of the gas may be 
ignored. Then with no pressure gradient required for the 
gas region, the pressure throughout the gas region has a 
uniform value. If the amplitude of the disturbance is small 
compared to its wavelength L, the pressure difference across 
the interface may be taken as approximately constant. Then 
the liquid pressure along the interface may be assumed to be 
approximately constant. 

The pressure gradient in the liquid approaches a 
constant value at some distance from the interface. 

However, at a distance up to a few wavelengths away from 
the interface, the pressure gradient is a function of x, 
and along a line y = constant has a larger magnitude in a 
region of nearest approach to the interface. Since the 
liquid velocity depends on the pressure gradient, the 
velocity of the liquid near the interface is greater in 
the region where the interface is leading than where it is 
lagging. Therefore, the portion of liquid which is lagging 
will continue to lag even further. Thus, the amplitude of 
the wave will continue to grow. Continued growth will 



result in long "fingers" of gas extending into the liquid. 

In the annular space between rotating cylinders, 
with an interface having a constant average position, the 
interface usually will oscillate axially about the mean 
position with a period equal to the time for one revolution 
of the rotating cylinder. This oscillation results from 
two unavoidable factors - eccentricity of mounting , and 
runout. Ignoring gravitational effects , the interface would 
be unstable only when moving toward the liquid. With the 
interface velocity in the opposite direction, the finite 
wave height, or finger length, would decrease. Therefore, 
fingers would be expected to form and grow to some finite 
length in a portion of the annulus , and then to become 
shorter and possibly disappear altogether in another portion 
of the annulus. 

The fingers thus formed would attain a maximum length 
which depends on the overall displacement of the oscillating 
interface . 

All of the fingers formed, as described above, would 
be expected to travel around the annulus at essentially 
half the surface speed of the rotating cylinder, although 
slight differences in their velocities could arise because 
of increasing or decreasing clearances along the path of 
fluid elements contained within such fingers, or variations 
in surface roughness, or other variations. Such velocity 
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differences could present an opportunity fpr a volume of 
air to become trapped within the liquid. If the top part 
of a liquid finger should travel with a tangential velocity 
slightly greater than that of the finger immediately preced- 
ing it, and for a sufficient length of time the gap could 
be closed near the top and thus trap a volume of gas beneath 
the interface. 

The chances of two adjacent fingers of liquid coming 
into contact with each other in this manner and trapping a 
volume of gas in the liquid would be enhanced by a high 
finger length-to -width ratio, i„e., by the presence of long, 
narrow fingers. Thus shorter wavelengths would be more 
conducive to trapping gas bubbles. Likewise, for a given 
wavelength, the greater the magnitude of the interface 
oscillation and thus the length to which a finger could 
grow, the greater the chance of entrainment of gas bubbles. 

The mathematical model shown in Figure 2, page 10, 
was established for an interface parallel to the x-axis , 

When the interface in this theoretical model has a velocity, 
the liquid is either advancing on both surfaces, or receding 
on both surfaces. However, if the physical model - a 
rotating cylinder inside a housing - has eccentricity but 
no runout, the interface at a particular angular position 
on the stationary surface would have a constant axial loca- 
tion. But the axial location would vary with angular 
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position on the stationary cylinder and therefore the 
interface would be oscillating axially relative to the 
rotating cylinder. The reverse would be true for runout 
of the rotating cylinder but concentricity of mounting. 
Therefore, it remains to be determined what conditions 
would produce an instability when an interface is advanc- 
ing or receding at different rates over the two surfaces 
between which the fluids are confined. 

A slight modification of the mathematical model 
shown in Figure 2, page 10, can be made to approximate the 
desired conditions. The plates in Figure 2 are shown to 
have equal but oppositely directed velocities U/2 in the 
x-direction. One plate is further assigned a velocity V 
in the y-direction, as shown in Figure 5. 

With this added condition imposed, appropriate 
boundary conditions are, 



at z = 0: 

u - -jfi 

u = + 

v = 0 ; 

(43) 


at z = c: 

v - V. 


The 

constants in 

equations 

(8) and (9) 

are unchanged, with 

the 

exception of 

C 2 , which 

becomes 



C . . V _ 
L 2 c 

1 (H + 

2y V 3y 

P g ) c., 

(44) 


and the y- component of velocity is 



Fluid 1 


cos Kx 


Fluid 2 


Figure 5. Model for analysis of a disturbed 
interface separating two viscous fluids moving transverse 
to the interface and confined between parallel plates 
having relative motion not parallel to the interface. 
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(45) 


Substituting equation (45) into equation (18) and 
integrating yields the average velocity, 


v = 


c 2 ft P 


12 y V 37 


fd P . \ , V. 

( 9y pg J 2 9y * 


(46) 


When equation (46) is integrated with respect to y, 
the solutions obtained for the pressures are 

p i * - (+i - w -$) - oi^ + c 7> t 47 ^ 

and 

^ 2p 2 / Vv 

P 2 - - — et (»2 - - i ) - ^ + C 8 * t 48 ) 

Equating the difference between the pressures from 
equations (47) and (48) with the pressure difference across 
the interface as given by equation (34) , noting that con- 
stant terms are eliminated as given in equation (39) , 
dividing by S and multiplying through by K = 2ir/L gives 
finally , 
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12m f -KS , 
— (y l e + 


w 2 e 


+ KS, 


2 7T 

L 




-y 2 ) + ( P]L - p 2 )g 


8ir 3 T 


(49) 


Equation (49) differs from equation (40) only in the 
coefficient of the viscosity difference. Other things 
being equal, if both the interface velocity and the plate 
velocity were directed toward the more viscous fluid, a 
larger magnitude of interface velocity would be required to 
cause the interface to be unstable to a small disturbance 
than would be required when neither plate moves vertically. 
If, for example, the vertical velocity of the plate and the 
interface velocity were equal and both moving toward the 
more viscous fluid, the value of interface velocity 
necessary to cause an instability would be twice the value 
required with neither plate moving vertically. Therefore, 
the physical model would have a greater tendency toward 
instability of the interface when both runout and eccen- 
tricity are present simultaneously. 

In assuming a pressure difference across the inter- 
face due to surface tension, to arrive at equation (34), it 
was assumed that P^ was greater than P 2 , and that the 
interface was concave upward. If the interface were concave 
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downward, the pressure difference would be given by 

CP 2 - Pi) - T(|- * (50) 

where R^ and R2 are radii of curvature having centers below 

the interface instead of above it, as assumed in obtaining 

equation (34). But the curvature, as approximated by 
2 7 

a S/3x , is positive when the center of curvature is above 
the interface. Therefore, in establishing a sign convention 
for a center of curvature located below the interface, the 
curvature is approximated by 


2 

1 3^S 



and equation (50) becomes 



or 




3 2 S \ 
3x 2 


(51) 


(52) 


(53) 


BK 2 e mt cos Kx 


The only difference between equation (54) and 
equation (37) is the algebraic sign on the constant term, 
which contains the clearance. But this term drops out 
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upon evaluation of the constants in the equations for the 
velocities. Therefore, the results are unchanged. 

For a gas -liquid interface in which the liquid is 
completely non-wetting, and the angle of contact is 180 
degrees, the conclusions concerning stability are identical 
to the conclusions reached for a liquid which wets the 
solid surfaces. 

III. STABILITY ANALYSIS - ACCELERATION EFFECTS 

The mathematical model for the analysis taking into 
account acceleration effects is shown in Figure 6. Viscous 
effects are neglected here and it is assumed that the fluid 
system undergoes an acceleration a normal to the interface 
and considered to be positive when directed downward, as 
shown in Figure 6. 

The coordinate system is assumed to move with the 
bulk fluid motion. Effectively, slip flow is assumed, and 
the distance separating the plates enters into the mathe- 
matical analysis only in the pressure difference across the 
interface, but does not appear in the final result. The 
initial mathematical formulation is essentially the same 
as that presented by Taylor [11] . 

The interface is assumed to have a small wave -like 
perturbation, with the equation of the surface of separation 
given by 



Figure 6. Model for analysis of the stability of 
an accelerated interface separating two fluids having 
different densities. 







32 


S = M- e mt cos Kx. (55) 

m 

The disturbance will be unstable when m is positive, and 
the mathematical procedure is similar to that followed in 
the viscous analysis. 

The instantaneous vertical component of the velocity 
of a small fluid element on the interface, relative to the 
mean position of the interface, is 


v = ff- = AK e mt cos Kx. (56) 

o T 

The appropriate velocity potentials are 

= - A e mt " Ky cos Kx, (57) 

and 

<j >2 = + A e mt + Ky cos Kx. (58) 


The equation of motion in the y- direct ion is 


9P = 
9y 


p(a - g 


DV^ 


9 V 
9 1 


- u 


.9 V 
9X 


v- 


(59) 


= p(a - g - 


9V 

9y 


(60) 
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The inertia terms are obtained from the velocity 
potentials in equation (57) by taking the appropriate 
partial derivatives, which are, 


u 


1 J F5T 


= KA e mt " Ky sin Kx , 


Vl . !ll . KA e mt - K y cos Kx, 


(61) 

(62) 


3x 


v 2. mt - Ky . 

- K A e ** sin Kx, 


(63) 


3V^ 

W~ 


- K 2 A e mt - Ky cos Kx, 


(64) 


3Vj 

FT 


= mKA e 


mt - Ky 


cos Kx. 


(65) 


Equation (59) becomes, for fluid 1, 


3P 

W 


I » p (a - g - mKA e mt * Ky cos Kx + Ku, 2 + Kv 2 ) . (66) 


For fluid 2, the required derivatives of the velocity 
potential are 


3 <J> 

u 2 “ 3x 


2 = - KA e mt + Ky sin Kx, 


(67) 
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^2 

w 


= KA e mt + Ky cos Kx, 


3v. 

i 

aT 


= - K 


2, mt + Ky . v 
A e 7 sm Kx, 


(6 8 ) 


(69) 



K 2 A e 


mt + 




8v 2 


at 


mKA e mt + Ky cos Kx * 


The equation o£ motion for fluid 2 becomes 
BP 

= P 2 ( a - g - mKA e mt + Ky cos Kx - Ku 2 2 - Kv 2 2 ) . 


Neglecting squares of the disturbance velocities 
equations (66) and (72) and integrating, the solutions 
the pressures are, 

P 1 = p i*^ a ' + P^A e mt " Ky cos Kx + C g , 


(70) 

(71) 


(72) 

in 

for 

(73) 




35 


P 7 = p ? (a - g)y - p 7 mA e 


mt + Ky 


cos Kx. 


2 - M 2 ¥ 2 

Assuming that the pressure difference across the 
interface, as given by equation (37), 


(74) 


t p i - p 2 ) - T (§* J) - 


(37) 


is applicable here, and setting equation (37) equal to the 
pressure difference obtained by subtracting equation (74) 
from equation (73) , gives 


(§- 


MI e” 1 cos Kx 
m 


) = ( - p l " p 2 ^ ~ g - ): 


+ mA e mt (p r e’ KS + p,e +KS ) cos Kx 


+ (C 9 - C 10 ) . (75) 

When cos Kx is zero, S is zero, .and equation (75) reduces 
to 

c“ T ^ C 9 ' C 1(P ‘ t - 76 - ) 

When equation (75) is divided through by S, as given by 
equation (55) , equation (75) becomes 

( Pl e _KS + P 2 e +KS ) = (a - g) (p 2 -p^ - TK 2 . (77) 



36 


2 

The algebraic sign of m is the same as the right- 

2 

hand side of equation (77). If m is negative, m is 
imaginary and the amplitude of the disturbance varies 
sinusoidally with time. This reflects the fact that 
viscous forces were neglected in the analysis . When m 
is positive , the sign of m could be either positive or 
negative » from a strictly mathematical standpoint. How- 
ever, the omission of Viscous forces in the analysis 
removes any means of damping, and it is therefore concluded 
from purely physical considerations that m will be positive 
when m is positive and the disturbance will increase ex- 
ponentially with time until, it attains an amplitude which 
is no longer small compared with the wavelength. 

When the gravitational acceleration acts toward the 
more dense fluid, equation (77) predicts that an instabil- 
ity will result from an acceleration of the fluid system 
toward the more dense fluid at a value which is sufficiently 
larger than gravity to cause the first term on the right to 
be larger than the secohd term. Therefore, surface tension 
contributes toward stability of the interface. 

The minimum wavelength of an unstable disturbance 
can be found from equation (77). The substitution K = 2 tt/L, 
where L is the wavelength of a disturbance, is made in 
equation (77) . 
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Unstable disturbances, having a wavelength L u> 
exist only when the right-hand side of equation (77) is 
greater than zero. Stated mathematically , 


4tt 2 T 


(a - g)(p 2 - Pl ) - > 0. 


( 78 ) 


u 


When equation (78) is solved for L u , the effect of surface 
tension is found to be to restrict uns table disturbances to 
those of wavelength greater than a critical value, given by 


u 


2 ttT 1 ^ 2 

(a - g) (p 2 


(79) 


Upon reaching a finite value, the amplitude of 


unstable disturbances will continue to grow and will result 
in the formation of fingers. The supporting reasoning is 
similar to that given for the growth of the unstable dis- 
turbances of the interface between the viscous fluids of 
Figure 4, page 21, into fingers. 

The model shown in Figure 4 is also applicable to an 
unstable interface in an accelerated fluid system. Consid- 
ering a gas-liquid interface) and ignoring the density of 
the gas, the interface pressure may be assumed constant. 
When the system is being accelerated toward the more dense 
fluid, the liquid in this case, the pressure in the liquid 
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at a location which is several wavelengths from the inter- 
face will be dependent only on the distance from the mean 
position of the interface. But at positions in the neigh- 
borhood of the interface, the pressure gradient near the 
leading portion of the interface will be greater than the 
pressure gradient near the lagging portion. Hence, the 
acceleration of the fluid near the leading portion will be 
greater in magnitude, and its velocity will become increas- 
ingly larger than the velocity of the fluid near the lagging 
portion. Thus, the amplitude of the disturbance will 
continue to grow. 

In an annulus between rotating cylinders, acceler- 
ations normal to the interface occur due to eccentricity 
and runout. Since the acceleration is alternately positive 
and negative, the wavelength of the fingers and the length 
to which fingers would grow would depend on the magnitude 
of eccentricity and runout, surface tension, shaft speed, 
etc., and would determine the potential for trapping gas 
bubbles in the liquid. 



CHAPTER III 


E XPE RI MENTAL STUDY 
I. EXPERIMENTAL FACILITY 

The test facility was designed and built to permit 
visual observations and photographic recording of an inter- 
face separating air and a liquid confined in the annular 
space between a rotating cylinder and a stationary housing. 

A transparent cylindrical housing is mounted to 
the test frame, as shown in Figure 7. Housings having 
nominal inside diameters of two, four, six or eight inches 
can be accommodated. Glass, and transparent acrylic 
plastic housings were used in the experimental work. An 
O-ring seal was used between the housing and the base. 

The stainless steel shaft, which is driven through 
a V-belt by a variable speed motor, was designed for mount- 
ing interchangeable test rotors ranging from two to eight 
inches in diameter. The test rotors, which were made of 
aluminum, have a length of six inches. The rotors are 
positioned on the shaft by a tapered shoulder on the shaft, 
and a tapered nut at the lower end. The rotors were ground 
and polished, and the runout was about 0.0003 to 0.0010 
inch . 
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Frame Rotating 
Mechanism — 



Frame 


Figure 7 . 


Interface study facility. 
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The shaft is supported by a ball bearing at the end 
on which the V-belt sheave is mounted. A helical groove 
hydrodynamic journal bearing having a diameter of one inch 
supports the other end, and is lubricated by the liquid 
being used in the interface study. 

Provision was made for rotating the test frame 
through an angle of 180 degrees, in order to attain any 
desired angle of orientation of the shaft axis in the 
gravitational field. 

With the frame orientation as shown in Figure 7, the 
volume below the test rotor, and within the housing, is the 
liquid reservoir. A hole was drilled and tapped in the 
base, and a valve was fitted into the passage to permit 
liquid to be introduced into the reservoir from below. 

This arrangement makes it possible to exclude all air 

t 

bubbles from the reservoir by filling it through this 
passage. The annulus is open to the atmosphere at the top 
of the rotor. 

A view of the test rig after completion of construc- 
tion is shown in Figure 8. The housing mounted in the test 
rig in this view is a five-inch outside diameter acrylic 
cylinder with an inside diameter of 4.116 inches. The 
aluminum rotor mounted on the shaft has a diameter of 
3.990 inches. In the left background of Figure 8 are other 
rotors having a nominal diameter of four inches, and also 
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Figure 8 


Experimental facilities 
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an extra shaft. The threaded rotors were used in viscoseal 
performance tests in a companion study. The white plastic 
funnel is attached to plastic tubing, which connects to the 
valve through which the reservoir is filled with liquid. 

In Figure 8, distilled water is in the reservoir, and the 
interface level is about midway up on the aluminum rotor . 

No provision was made for cooling the working fluid. 
However, a thermocouple was installed in a radial hole 
drilled through the acrylic cylinder to provide for moni- 
toring the temperature of the liquid in the annulus. 

Visual observations were facilitated by the use of 
a stroboscope, which was also used as the light source for 
still photographs of about 1.2 microseconds exposure time. 
The still photographs were made with a 35 mm camera using 
black and white film. 

High speed movies were made at film speeds up to 
5,000 frames per second using black and white film. 
Illumination was provided by two lamps, each having a 
rating of 1,000 watts. 

II. DESCRIPTION OF EXPERIMENTAL INVESTIGATION 

An annulus with a nominal diameter of four inches 
was arbitrarily chosen for the initial experimental studies, 
which preceded the analytical development presented in 
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Chapter II. Because of the nature of the subsequent ana- 
lytical development, it was not found necessary to employ 
the full range of diameters for which the facility was 
designed. Consequently, annuli with a nominal diameter of 
four inches, but with a variation in the thickness, were 
employed throughout the experimental study . 

Three radial clearances were used in the experimental 
work. Specific dimensions of the rotor-housing combinations 
are given in Table I. 


TABLE I. 

DIMENSIONS OF ROTOR-HOUSING COMBINATIONS 
USED IN THE EXPERIMENTAL WORK 


Rotor Diameter 
(Inches) 

Housing I.D. 
(Inches) 

Avg. Radial 
Clearance (Inches) 

3.990 

4.002 

0.006 

3.980 

4.002 

0.011 

3.990 

4.116 

0.063 


A 0.0 06 -inch average radial clearance was provided by 
two different housings, one made of glass and one made of 
acrylic, and each having an inside diameter of 4.002 inches. 

In the early experimental work, which was of an ex- 
ploratory nature, air bubbles were observed to be present 
in the annulus at sufficiently high shaft speeds. An effort 
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was made to establish an empirical criterion for the onset 
of air entrainment. Following the approach taken by McGrew 
and McHugh [5], and by King [6], the Reynolds number and the 
Weber number were considered in attempting to obtain a cor- 
relation with the beginning of air entrainment. 

Studies were made at shaft speeds up to 4,000 rpm 
using fluids having viscosities from 0.5 centipoise for 
methyl alcohol up to 1200 centipoise for glycerine. Values 
of surface tension for the fluids employed ranged from 22 
dynes/cm. for ethyl alcohol to 72 dynes/cm. for water. No 
significant relation was established between the beginning 
of air entrainment and the Weber number or Reynolds number, 
or combination of the two numbers. However, it was found 
that air entrainment occurred for some fluids, particularly 
the more viscous fluids, under laminar flow conditions. For 
a Reynolds number based on the surface velocity of the rotor 
and on the radial clearance, Couette flow remains laminar 
for Reynolds numbers below 1,000 [14] , and may remain laminar 
at Reynolds numbers of four times this value [15] . Air en- 
trainment was observed to occur at Reynolds numbers of 100 
or less, placing the flow well into the laminar regime. 

The observation that air entrainment could occur under 
laminar flow conditions, as well as under turbulent flow 
conditions, provided insight as to the general direction 



46 


which the analytical study should take. This eventually led 
to the theoretical prediction of an interface instability 
not closely related to either the Reynolds number or the 
Weber number. 

Experimental verification of the predicted instabil- 
ity of the interface was primarily of a qualitative nature 
only. While the runout and eccentricity in the experimental 
facility provided the conditions to produce an interface in- 
stability, they also presented problems in determining local 
velocities and radial clearances. 

However, visual observation revealed an oscillating 
interface, and by varying the shaft speed and the viscosity 
of the liquid used, conditions could be produced which would 
cause the formation of fingers of air extending into the 
liquid. These fingers were seen with the aid of a strobo- 
scope, and air bubbles trapped below the interface were also 
seen. But the actual development of the fingers and the 
trapping of the air bubbles were not seen satisfactorily, 
except at very low shaft speeds . 

For the higher shaft speeds, a high-speed movie 
camera was used to facilitate viewing the development of 
fingers and the formation of air bubbles. Film speeds up to 
5,000 frames per second were used for shaft speeds up to 
2140 rpm. 



CHAPTER IV 


EXPERIMENTAL RESULTS 

Air entrainment occurred with each of the three 
clearances employed in the experimental study. The shaft 
speed at which air entrainment occurred varied with the 
thickness of the annulus, and with the properties of the 
liquid being used. 

Clearance of 0.006 inch 

The results obtained with a glass housing providing 
a clearance of 0.006 inch were not identical to the results 
obtained with an acrylic housing providing the same clear- 
ance, but the difference in the results was not attributed to 
the difference in the materials of which the housings were 
made. There were differences of comparable magnitude between 
different test runs using the same acrylic housing. The dif- 
ferences were attributed primarily to unavoidable variations 
in runout and eccentricity. These variations were created 
when the apparatus was disassembled for cleaning, or other 
purposes, and then reassembled with a slightly different 
alignment of rotor and housing. 

The eccentricity ratio, which could not be precisely 
determined, was estimated to be about 0.1 to 0.4 for tests 
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conducted with the glass housing and with the acrylic 
housing . 

Using a 50% glycerine and 50% water solution, having 
a viscosity of five centipoise at room temperature, the 
formation of pronounced fingers was observed at relatively 
low shaft speeds. The photograph shown in Figure 9 was 
taken when the shaft speed was 50 rpm. The corresponding 
Reynolds number was eight, indicating the presence Of 
laminar flow conditions. This photograph depicts a typical 
interface profile with pronounced fingers formed under the 
influence of viscous effects. 

Using the same solution, a few hundred feet of 16 
millimeter movie film were taken for shaft speeds up to 
600 rpm. The progressive stages of the development of the 
fingers can be seen on this film. A typical movie film 
sequence is shown in Figure 10 for a shaft speed of 440 rpm. 

Interface details were not easily discernible when 
viewing the motion picture projected at normal speeds, and 
were quite difficult to see in single frame viewing. For 
this reason, the same film sequence is reproduced as Figure 
11, where the interface has been traced with ink on the 
photographic prints in order to assist in visually following 
the changes in the interface seen in the unaltered photo- 
graph of Figure 10. 
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Figure 9. Fingers of air and 50% mixture of 
glycerine and water at shaft speed of 50 rpm. 



E 







Figure 10. Motion picture frame sequence showing 
development of fingers in 50% glycerine-50% water solution 
Shaft rotation of 26.5° between frames, shaft speed - 440 
rpm, clearance - 0.006 inch. 



Figure 11. Retouched motion picture frame sequence 
showing development of fingers in 50% glycerine -50% water 
solution. Shaft rotation of 26.5° between frames, shaft 
speed - 440 rpm. clearance - 0.006 inches. 
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The film speed for the movie sequence shown was 500 
frames per second, and every fifth frame of the film is shown 
in Figures 10 and 11. Elapsed time between adjacent frames 
shown was 0.01 second, representing the time for 26.5 degrees 
of shaft revolution. 

The interface was descending in the first nine or ten 
frames shown. This is not immediately apparent in the photo- 
graphs shown, but was readily seen in viewing the projection 
of the film onto a screen at normal projection speeds. 

A dark vertical line and a dark horizontal line, 
each passing through the center of the picture, are seen 
in each frame. These reference lines were built into the 
optical system of the camera. To the right of and adjacent 
to the vertical reference line is seen a light vertical 
band, which is one of the bolts holding the plate which 
clamps the transparent housing into position. The bolt was 
wrapped with white paper, and a film identification number 
was taped to it. 

A dark band is seen to the left of the vertical 
reference line. This is a section of the aluminum rotor 
which was insufficiently lighted because of the difficulty 
in uniformly illuminating a curved, highly reflective sur- 
face. 

When using distilled water in the annulus, severe 
gas ingestion generally occurred at shaft speeds above 
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2150 rpm , and interface irregularities were observed at 
shaft speeds above 1100 to 1350 rpm. The variation in 
minimum speed at which these irregularities became evident 
was attributed to the unavoidable differences in eccentric- 
ity and runout previously discussed. 

On one occasion, with the shaft operating at 1100 
rpm, the interface showed no irregularities and no trapped 
air bubbles were visible. As the shaft continued to operate 
at 1100 rpm, the shaft axis was rotated 180 degrees, placing 
the liquid above the interface and the air below it. The 
shaft continued to operate at this speed, with the interface 
inverted, for over an hour with no visible evidence of air 
bubbles in the liquid or any loss of liquid from the inter- 
face . 

However, on another occasion, when the shaft eccen- 
tricity was probably different, interruptions of the 
interface and air ingestion occurred for the same shaft 
speed with the shaft oriented in the normal position (inter- 
face not inverted). This is shown in Figure 12. 

A comparison of Figure 12 with Figure 9, page 49, 
indicates interface irregularities and trapped air bubbles 
of smaller dimensions at the higher shaft speed. Although 
viscous effects and acceleration effects on the interface 
were present simultaneously, the viscosity of the water in 
Figure 12 was only 20 per cent of the viscosity of the 



Figure 12. Air-distilled water interface region 
at shaft speed of 1100 rpm. 





glycerine -water solution in Figure 9 and suggests a stronger 
influence of acceleration effects on the air entrainment in 
Figure 12. 

Clearance of 0.011 inch 

The eccentricity ratio was estimated to be less than 
0.2 for observations made with a clearance of 0 .011 inch. 

Using water in the annulus, air bubbles were present 
in the annulus for shaft speeds above 1100 rpm. 

The shaft axis was rotated 180 degrees while the 
shaft was operating, and no leakage from the interface was 
ob s e rve d . 

Clearance of 0.063 inch 

An eccentricity ratio of less than 0.1 was estimated 
for a clearance of 0.063 inch. 

Using water in the annulus, air bubbles appeared in 
the annulus at shaft speeds of 670 rpm or higher. The 
bubbles were observed to be aligned in circumferential rows 
with an axial spacing varying from about 0.15 inch at low 
speeds up to about 0,25 inch at speeds above 3000 rpm. 

The rows of bubbles had not been anticipated, and 
aluminum powder was mixed with the water in the annulus in 
order to verify that the spacing of the rows of bubbles 
corresponded to the spacing of vortex pairs predicted by 
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Taylor [16]. The initiation, theoretically and experimen- 
tally, of the Taylor instability was found to be 25 rpm. 

With increased shaft speed, the axial spacing of the 
Taylor's vortex pairs was found to increase from the 
0.125 -inch spacing predicted by Taylor for the onset of 
the instability. 

Efforts to operate the shaft in the inverted position 
without the loss of water were unsuccessful. A large range 
of shaft speeds was employed in attempts to rotate the shaft 
axis, while the shaft was rotating, without a breakup of 
the interface. In all efforts, water was lost from the 
annulus before the shaft axis was rotated through an angle 
of more than about 120 degrees. 



CHAPTER V 
DISCUSSION 

The theoretical development presented in Chapter II 
and experimentally verified in a qualitative way, with 
typical results displayed in Chapter IV, provides an expla- 
nation of a possible mechanism by which gas ingestion can 
occur in the region of a gas-liquid interface. 

The operation of a viscoseal involves regions in 
which the interface has a component of velocity relative to 
one of the surfaces such that the liquid is advancing on a 
solid surface, and other regions where the interface is 
receding on the solid surface. With runout of the rotating 
shaft, there may be regions in which the liquid is receding 
on both surfaces. An interface which is receding on one or 
both surfaces may be subject to the viscous type of inter- 
face instability, which could lead to gas entrainment. 

Fluid particles adjacent to the interface in the 
viscoseal are also subjected to local accelerations having 
a component normal to the interface* These normal acceler- 
ation components, when directed toward the liquid, can also 
cause an interface instability resulting in gas entrainment. 

The subsequent motion of gas bubbles thus trapped 
would depend on a number of variables, including the orienta- 
tion of the viscoseal in the gravitational field, the 



viscoseal geometry, whether turbulence exists, and other 
factors. If all entrained bubbles subsequently moved back 
to the interface rather than being transported away from the 
interface, the effect of the gas entrainment could be incon- 
sequential with respect to the satisfactory operation of the 
viscoseal. However, if entrained gas bubbles continually 
accumulate in the liquid confined in the annular space, the 
viscoseal would most frequently be rendered ineffective to 
its intended purpose. As. discussed in Chapter I, the 
mechanics of the gas transport process is a separate aspect 
of the problem associated with gas ingestion, and is not 
well understood. 

In the viscoseal, a small element of fluid adjacent 
to the interface has motion with respect to the grooved 
surface as well as having motion with respect to the smooth 
surface. An unstable disturbance, in the form of a 
developing finger, will have the same general motion as 
the local average motion of the liquid element in which 
the disturbance is located. An undetermined but finite 
amount of time is required for a small disturbance to grow 
to a sufficient size to result in the joining of the tips 
of two liquid fingers to isolate a quantity of gas as a 
bubble. If, during the time of growth of such a disturbance, 
the liquid element and its developing disturbance move into 
a region which is stable to a small perturbation, then 
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the disturbance will no longer continue to grow, but will 
begin to decrease in size and no gas bubble will have 
resulted from the disturbance. 

Therefore, the theory presented here does not lead to 
a quantitative prediction of the actual trapping of a gas 
bubble in the liquid. Thus regimes of viscoseal operation 
in which gas ingestion would be expected could not be deter- 
mined, even if the motion of the fluid field were known. 

However, if the fluid motion were known, operating 
ranges which could be definitely expected to be free of gas 
ingestion could be established by determining that interface 
stability existed on all parts of the interface during the 
entire time of one shaft revolution. This assumes that 
there are no mechanisms of gas ingestion other than those 
presented here. 

Other physical applications approximately fitting the 
hypothetical model include face seals, buffered bushing 
seals, and journal bearings. 

Journal bearings are inherently vulnerable to gas 
entrainment resulting from the type of interface instability 
presented here, since any radial load results in eccentric 
operation of the journal. If the bearing is not immersed in 
the lubricant, and a gas-liquid interface exists in the 
annulus, then over a portion of the interface the liquid is 
receding from the solid surface of the shaft. This 
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application represents rather closely the model on which 
the analytical study was based, and the journal bearing 
operation is essentially identical to the experimental 
facility used in the interface stability study. However, 
cavitation (where the local pressure becomes less than the 
vapor pressure of the lubricant) could occur in the annulus 
of a journal bearing which is completely immersed and in- 
volves a different mechanism than that causing what is 
referred to here as gas ingestion. 

Buffered bushing seals are similar in many respects 
to the journal bearing, and gas ingestion encountered in the 
buffered bushing seal may well have its origin in the same 
type of interface instability which is postulated here to 
occur in the journal bearing. 

Face seals represent, ideally, relative motion between 
two closely spaced parallel plates between which a gas-liquid 
interface exists. Due to non -uniformity of the spacing be- 
tween the plates, the interface would not have the shape of 
a circle and the distance from the interface to the center 
of rotation would vary along the interface. Therefore, a 
receding interface would exist in a region or regions, and 
fluid elements in the interface would also undergo normal 
accelerations. Hence, potential conditions for instability 
of the interface and for gas ingestion are present. The 
potential for gas ingestion would appear to be considerably 
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greater when sealing a liquid in the center from a gas, 
since the centrifugal acceleration of the interface is 
directed toward the liquid. 

Helically-grooved face seals have characteristics 
similar to the viscoseal , in that the fluid flow in the 
vicinity of the interface is considerably complicated by 
the presence of the grooves and lands. Those designed for 
inward pumping of the liquid would appear to have a greater 
tendency toward gas ingestion than face seals with plain 
surfaces . 

The results in the experimental work were frequently 
not satisfactorily reproducible. As discussed earlier, some 
of the variations in results were attributed to differences 
in shaft and housing alignment. However, there were also 
occasions when there was no reason to believe that there 
were changes in dimensions of the physical equipment, and 
yet the gas ingestion characteristics varied substantially 
from one observation to another. These variations were 
particularly noticeable when there was a time interval of 
several hours between an observed event and an attempt to 
duplicate the event. 

The mathematical model on which the theoretical 
development was based was considerably simplified in order 
to permit an analytical solution. A contact angle of zero 
was assumed in order to determine an expression for the 
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pressure difference across the interface. It was noted 
in Chapter II that the mathematical result is unchanged by 
assuming a non-wetting interface, and it is seen that for 
any uniformly constant interface contact angle the mathemat- 
ical results are unchanged. 

However, Dettre and Johnson [13] have shown that 
advancing and receding contact angles may differ by as much 
as 80 degrees. Also, nonuniformity of the surface causes 
local variations in contact angles. These authors also 
imply that the relative humidity of the air in contact with 
an air-liquid interface affects the angle of contact, and the 
time of exposure of a liquid to a surface is stated to in- 
fluence the contact angle. 

In the experimental work reported here, the relative 
humidity of the air in the laboratory was not controlled, 
nor was it recorded. Also, changes in the physical appear- 
ance of the surface of the aluminum rotors were observed, 
causing speculation that oxidation or other changes were 
producing changes in the wetting characteristics on the sur- 
face . 

One particular variation in properties observed to 
result from long time intervals between experimental obser- 
vations was the composition of mixtures of glycerine and 
water, with the attendant variations in viscosity and other 
properties. It was discovered that the amount of water in 
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solution increases or decreases significantly when the 
solution is not in equilibrium with the air to which the 
liquid surface is exposed. Then after equilibrium is 
reached, measurable changes in composition occur from one 
day to the next as the relative humidity of the air increases 
or decreases. The room temperature was held within a narrow 
range, winter and summer, but the relative humidity was not . 

These factors are pointed out to indicate some of the 
variables which may enter into the problem of obtaining 
reproducible results in viscoseal operation, or other related 
applications . 

Although the problem of reproducibility was 
experienced in this work, the problem did not prevent the 
achievement of the objective of the tests - that of verify- 
ing the mechanism by which gas ingestion can occur in an 
interface between rotating cylinders. 



CHAPTER VI 


CONCLUSIONS AND RECOMMENDATIONS 

Some conclusions evolving from this investigation 
are that : 

1. A phenomenological explanation of gas ingestion 
can be obtained through a stability analysis of a dynamic 
gas -liquid interface, 

2. For fluids contained between closely spaced 
surfaces, gas ingestion can be caused by a gas -liquid 
interface instability resulting from a velocity and/or an 
acceleration of the interface toward the liquid. 

3. Surface tension tends to prevent interface 
instabilities and gas ingestion. 

4. Gas ingestion in viscoseals appears to be an 
essentially inherent phenomenon which can not be readily 
eliminated through optimization of the geometry. 

5. Eccentricity, shaft runout and vibrations con- 
tribute to conditions necessary for the instability of an 
interface in rotating machinery. 

Although a mechanism of gas ingestion has been 
established, both theoretically and experimentally, this 
study did not provide a quantitative criterion for predict- 
ing gas ingestion. For specific engineering applications 
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in which a gas-liquid interface must be present in an 
annulus between rotating shafts, information is needed to 
aid in selecting design parameters which will provide 
operation free of the problem of gas ingestion. 

Since gas ingestion in viscoseals seems to be a 
phenomenon which is not easily avoided, efforts should be 
made to design a method of separating out the ingested gas. 
A substantial effort in this direction has been reported 
by Boon [17]. A further effort to provide a recirculation 
and centrifugal separation path in the viscoseal is being 
made by the Mechanical and Aerospace Engineering Department 
at The University of Tennessee. 
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